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A continuous sequence of infinitesimal unitary transformations is used to diagonalize the 
quantum sine-Gordon model for G (27r,cx3). This approach can be understood as an exten- 
. sion of perturbative scaling theory since it links weak- to strong-coupling behavior in a systematic 

expansion: a small expansion parameter is identified and this parameter remains small throughout 
the entire fiow unlike the diverging running coupling constant of perturbative scaling. Our approxi- 
mation consists in neglecting higher orders in this small parameter. We find very accurate results for 
the single-particle/hole spectrum in the strong-coupling phase and can describe the full crossover 
from weak to strong-coupling. The integrable structure of the sine-Gordon model is not used in our 
approach. Our new method should be of interest for the investigation of nonintegrable perturbations 
and for other strong-coupling problems. 
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B] 

• I. INTRODUCTION 

> 

A. Motivation 

■ 

^ . Perturbative scaling theory plays a key role for analyzing the large class of physical systems with a mismatch between 
^ ■ the high-energy scale of the model and the experimentally interesting low-energy scale. For example, in field theory 
(— ( one is generally interested in the universal properties at energies much lower than the UV (ultraviolet)-cutoff, or 
Q in condensed matter physics in energies and temperatures much smaller than the Fermi energy/temperature usually 
O of order of a few 1000 K. In order to link high-energy and low-energy regimes it is of fundamental importance 
to perform perturbation theory in a stable order by first analyzing the effect from high-energy scales, and then 
I ' progressively smaller energies. An elementary example for this procedure is provided by atomic physics, where one 
^ , e.g. first establishes the fine structure of a spectrum before using these states to evaluate the hyperfine splittings. 
CO ' For systems with continuous energy scales, like in field theory, the above observations have led to the development 
' of perturbative scaling theory. Perturbative scaling ideas have become a key theoretical tool for analyzing physical 
' systems with many degrees of freedom. The principal idea is to study perturbatively the effect of lowering the high- 
energy cutoff by finding a Hamiltonian with this reduced cutoff and renormalized couplings that describe the same 
low-energy physics as the original Hamiltonian. In a path integral formulation this is conveniently achieved by 
■ successively integrating out the high-energy degrees of freedom, 
-j—; This procedure leads to the well-known renormalization group (RG) equations that describe the flow of the running 

coupling constants upon lowering the UV-cutoff. For the important class of strong-coupling problems, however, 
^ the RG-equations lead to running coupling constants that grow larger and larger at smaller energy scales (and 
I often eventually even diverge). Since the RG-equations themselves are derived perturbatively, this means that the 
■ O perturbative scaling approach breaks down for strong-coupling problems. Well-known examples for this class of 
^ ■ models are the Kondo model in condensed matter physics or QCD in elementary particle physics. In spite of its 
eventual breakdown, perturbative scaling can still contribute significantly to the understanding of strong-coupling 
problems. For example, in the Kondo model, the divergence of the running coupling constant occurs at an energy 
. !^ thai sets the low-energy Kondo scale of the model, which already allows considerable insight into the problem. Still 
K> , the approach becomes uncontrolled since the coupling constants grow very large, and it has, so far, not been possible 
5^ ' to extend the perturbative scaling approach in such a way that a controlled systematic expansion emerges that links 
I weak- to strong-coupling behavior. One can sum up these observations by noting that the perturbative scaling 
approach often allows us to identify the relevant low-energy scale of a strong-coupling problem, but frequently not 
the physical behavior associated with this energy scale or the crossover behavior linking high and low energies. For 
an excellent review of these issues see Ref. [|] . 



*This work is dedicated to Prof. Wegner on the occasion of his 60th birthday. 
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This paper will exemplify the way in which Wegner's method of flow equations can overcome these shortcomings 
and provide an analytic description for a weak- to strong-coupling behavior crossover. In the flow equation approach, 
a continuous sequence of infinitesimal unitary transformations is applied to a many-particle Hamiltonian such that 
the Hamiltonian becomes successively more diagonal. Wegner has set up this approach in a differential formulation 



Here ri{B) — —rj{B)^ is an anti-Hermitian operator. Therefore H{B) as obtained by the solution of this differential 
equation describes a one-parameter family of unitarily equivalent Hamiltonians. H{B = 0) = iJ is the initial condition 
relating us to the original Hamiltonian H in which we are interested. We want H{B — oo) to he diagonal. In order 
to achieve this, Wegner has proposed a suitable choice for the generator rj{B) that we will discuss in more detail 
in Sect. HI. A. Wegner's construction of 77 generates a Hamiltonian flow where the interaction matrix elements that 
couple degrees of freedoms with a large energy difference are removed first (for smaller B), and more degenerate matrix 
elements during later stages of the flow. This separation of energy scales is reminiscent of the perturbative scaling 
approach and allows a stable sequence of approximations. As opposed to the perturbative scaling approach, however, 
degrees of freedom are not integrated out in the flow equation approach, instead they are successively diagonalized. 
A similar framework that contains Wegner's flow equations as a special case has independently been developed by 
Glazek and Wilson {similarity renormalization scheme) 

So far the flow equation approach has been applied to a variety of models in condensed matter theory like the 
n-orbital model 0,0] , impurity models like the spin-boson model g and the Anderson impurity model |^ , electron- 
phonon systems and spin models |^,|lO| etc. (for an overview see also Ref. ^^). One advantage of this scheme 
lies in the observation that it is a non-perturbative approach due to the separation of energy scales, but still has 
access to all energy scales since no degrees of freedom are integrated out. Therefore one can investigate correlation 
functions on all energy scales Q. Also the flow equation approach allows the systematic derivation of low-energy 
effective Hamiltonians not plagued with singular interactions that frequently occur in other approaches 

However, these applications did not deal with strong-coupling problems as defined above, which would be a very 
interesting perspective for this new method. Glazek and Wilson undertook a first step in this direction in Ref. [Q. 
They investigated a quadratic Hamiltonian that shows strong-coupling behavior due to the formation of a bound 
state from a continuum, and demonstrated how this model can be solved using infinitesimal unitary transformations. 
However, since they dealt with a quadratic Hamiltonian, this was not a true many-particle strong-coupling problem 
as would be of most interest in condensed matter theory or high energy theory. 

Recently, I described the application of the fiow equation method to the one-dimensional quantum sine-Gordon 
model [ p^ . The sine-Gordon model is a many-particle problem with an interesting phase structure including a 
strong-coupling regime. It was shown in Ref. jT^ that it is possible to use the fiow equation scheme to develop a 
systematic expansion that links weak- to strong-coupling behavior in a controlled way. Already the leading order of 
this expansion was is close agreement with exact results. In the present paper I will present the various details of the 
calculation not included in the original Letter in a self-contained manner. 

The sine-Gordon model is defined by the Hamiltonian 



where (/'(x) is a bosonic field and n(x) its conjugate momentum field with the commutator [H(a;), 4>{x')\ = —iS{x — x'). 
M > is a small dimensionless coupling constant and A cx r^^ an implicit UV-cutoff. We are interested in the universal 
properties for energies \E\ <^ A. 

The sine-Gordon model exhibits a strong-coupling phase for ^ Stt with a mass gap and fermionic low-energy 
excitations {massive solitons). The perturbative scaling analysis leads to a characteristic strong-coupling divergence 
of the running coupling u in this regime. This makes the sine-Gordon model an interesting test model for our new 
approach. The main emphasis in this paper will not lie in deriving new results in this well-studied model, but in 
showing how these results follow within the flow equation method, and how therefore our new method can be useful 
for strong-coupling problems more generally. 

Other features that make the sine-Gordon model an attractive test model are its interesting phase structure with 
a Kosterlitz-Thouless type transition to a phase with massless solitons at /3^/87r « 1 + 0{u), its integrable structure 
that allows the comparison with exact results and its relation to a variety of other models like the spin-1/2 

X-Y-Z chain, the Id electron gas with backward scattering, the Thirring model in held theory and the 2d Coulomb 
gas (for an overview of these relations see Ref. jl^). Therefore the results from the flow equation approach can be 
viewed within a variety of model contexts. 



^[rj{B),H{B)] . 



(1) 




(2) 
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The main motivation for being interested in the flow equation approach to this integrable model hes, however, 
in the observation that our new method does not make use of the integrable structure. In our approach a small 
parameter is identified and used within a suitably renormalized perturbation expansion. The usual perturbative 
scaling approach fails because the initially small expansion parameter u diverges during the RG-procedure. In the 
flow equation approach the expansion parameter will turn out to be the product of the running coupling u and a factor 
(— 1 + /3^/47r). This combination will always remain small during the flow. It is therefore feasible to study for example 
nonintegrable perturbations and correlation functions within our new approach, which should be of considerable 
interest in a variety of contexts. Although the calculations presented here appear rather lengthy and technical at 
first, they are straightforward and much closer to conventional many-body techniques than methods building on the 
integrable structure. 

B. Outline 

The structure of this paper is as follows. Sect. II deals with some general properties of the sine-Gordon model 
that are important in the sequel. In Sect. II. A the sine-Gordon model and the regularization used in this paper 
are introduced. Sect. II. B reviews the perturbative scaling analysis, the phase structure, and the strong-coupling 
behavior. In Sect. II. C various exact results based on the integrable structure of the sine-Gordon model are summed 
up, especially properties of the point — Air where the model becomes equivalent to a noninteracting Thirring model. 
This equivalence will play an important role in understanding the structure of our flow equation approach later on. 

After setting the stage in Sect. II, Sect. Ill deals with the actual application of the flow equation approach to the 
sine-Gordon model. Some general properties of the flow equation method are reviewed in III. A. Then the appropriate 
generator r]{B) for the sine-Gordon model diagonalization is worked out in III.B and the commutator [ri{B),Ho] 
evaluated in III.C. The key computational parts of the flow equation approach are contained in III.D and III.E, where 
the commutators [ri{B), Hiat{B)] and [ri{B), Hdia.g{B)] are evaluated. From these commutators the flow of 0^{B) and 
of the running coupling constant are deduced in Sect. III.F. 

For the convenience of the reader, all the results from this technical part are summed up in Sect. IV. A, in particular 
the Hamiltonian H{B) along the flow and the set of flow equations governing the various parameters in H{B). We 
will see in IV. B that in the strong-coupling phase H{B) flows to an effective low-energy noninteracting Thirring 
model. The mass gap of the sine-Gordon model can be easily deduced from this low-energy model, and the results 
are then compared with perturbative scaling analysis and exact integrable model results. The agreement will turn out 
to be very good. In Sect. IV. C the final diagonal Hamiltonian H{B = oo) is discussed in more detail, in particular 
the soliton dispersion relation and properties in the crossover region. Finally in Sect. IV. D the approximations and 
the expansion parameter of our approach are reviewed. 

Sect. V sums up the conclusions and an outlook to open questions. The Appendix contains important properties 
of vertex operators that are used throughout this paper. 

II. SINE-GORDON MODEL 
A. Definition 

The one-dimensional quantum sine-Gordon model is defined by the Hamiltonian 

H = Jdx + i (^1^^ cosmx)]^ . (3) 

4>{x) is a bosonic field and n(a;) its conjugate momentum field with the fundamental commutator 

[U{x),(l){x')]^-id{x-x') . (4) 

In an UV-momentum cutoff A oc t"-^ is implied, u is a dimensionless coupling constant. Without loss of generality 
we will assume u > and /3 > 0. 

Expanding the fields in normal modes gives 
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-^E#e"^'^(^iW+'^2(fc)) (5) 



Sums over wavevectors k,p,q, . . . are to be understood in the sense 



n(x) = v4^e-'=^ (ai(fc) - . (6) 



fc n= — oo 

with fc = 2Tm/ L throughout this paper. L is the system size. The basic commutators for k,k' > are 

[ai(-fc),ai(fc')] - ['T2(fc),fT2(-fc')] -4fe' 

[a,(fc),a,(fc')] = K-(-fc),<7,(-fc')]-0 (8) 

and for j ^ j' 

[a,(fc),a,v(fc')] = h(-fc),a,,(fc')] = [a,(-fc),c7,,(-fc')] = . (9) 
An other commutators can be derived via (t|(— fc) = (k). The vacuum \n) is defined by 

ai{-k)\n) = a2ik)\n) = (10) 
for all fc > 0. The notion of the dual field 9(x) will also be useful. Q{x) is defined by 

a,e(x) = -n(a:) , (11) 

leading to the commutator 

[Q{x),(l){x')]^ie{x-x') . (12) 

In terms of normal modes one finds 

Six) = --^E 4^ (^i(fc) - '^^(fc)) • (13) 

The concept of vertex operators will play an important role in the sequel. Vertex operators Vj (a; a;) are defined as 
normal-ordered exponentials 



Vj (a; x) = : exp j ±a V^M g 



1^1 ^—^\p\-'ipx 



a,ip) : (14) 



with + (upper sign) corresponding to j = 1 and — (lower sign) to j = 2. This sign convention will be used throughout 
this paper. Normal ordering : . . . : amounts to commuting all the operators that annihilate the vacuum according to 
(O) to the right. One can rewrite dlj) in terms of the field and its dual (|ll|) 



Vj{a; x) — : exp ^iay/ir J dec{e) [(f>{x + e) ± 8(a; + e)]^ : 



(15) 



with the Lorentzian 



. . a/2TT 

+ 0-^/4 



c(e) is normalized 

POO 

dec{e) = l (17) 
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and c(e) '^—^ S(e). Further properties of vertex operators, in particular their operator product expansion (OPE), are 
reviewed in the Appendix. 

One can rewrite the interaction term of the sinc-Gordon model (||) in terms of vertex operators 



H 



dx 



-tf(.x) 



1 

2\dx 



27ra^ 



[Viia- x)V2{~a] x) + V2{a; x)Vi (-a; x)) 



(18) 



The prefactor {2na/L)°' follows from ([Aq). Here and in the sequel a and /3 are used interchangeably with the 
identification 



dcf (3 



An 



(19) 



Eq. (18) is the form of the sine-Gordon Hamiltonian that we will investigate with the flow equation approach: No 
implicit momentum cutoff is implied in dlj) : The UV-regularization of the Hamiltonian (|l^) is achieved by the cutoff 



parameter a > in the vertex operators. The regularizations in ( |18| ) and (g) are related by a ocAcxr . For a 
direct comparison between the original Hamiltonian (||) and the form (|l^) used here one can also identically express 
(H as 



H 



dx 



-n^x) 



dx 



[3 I dec{e)(f){x + e) 



(20) 



The regularization with a therefore amounts to smearing out the interaction term. In the limit a — > one recovers 
the cos(/3</)(2;))-interaction term. 

The universal properties of the sine-Gordon model for energies \E\ <C a^^ are not affected by this choice of 
regularization. We find, however, notational simplifications and more compact expressions in the course of our 
calculation when we start with (^8|) (or equivalently (pO|)). In order to clarify the main conceptual ideas of the flow 
equation approach, it will therefore be convenient for us to use the regularization ( |T8| ) with the UV-cutoff a^^ built 
in via the definition of the vertex operators. 



B. Perturbative scaling analysis 

The flow equation approach can be viewed as an extension of perturbative scaling. Therefore it is useful to briefly 
review the results of the perturbative scaling analysis as applied to the sine-Gordon model. A comprehensive review 
can be found in Ref. JlTt . 

In 2-loop order there are two renormalization group equations that describe the flow of u and /3 upon integrating 



out the degrees of freedom with A — dA < |fc| < A (see Ref. |18 

dp-^ _ 
dlnA ~ 47r 



''1^-21.. (21) 



dhiA yi-K 

The initial conditions are u(r^^) = u and /3(t^^) = /3. These scaling equations give rise to the Kosterlitz-Thouless 
type phase diagram shown in Fig. ^ The two separatrices S± originating from — Sir with (3"^ — 87r(l ± u) for 
small u divide the parameter space in three sectors: 

1. the weak-coupling sector I; 

2. the crossover sector II; 

3. the asymptotic freedom sector HI. 

Both in II and III the perturbative scaling equations ( pT| ) lead to strong-coupling behavior with the running coupling 
constant u growing larger and larger during the flow. Therefore the perturbative RG approach eventually becomes 
invalid in these sectors. This indicates the opening of mass gap in the spectrum in II and HI. Across S+ the system 
undergoes a Kosterlitz-Thouless type phase transition between this massive phase and the massless phase in I. 



5 




FIG. 1. Perturbative scaling flow and the phase diagram of the sine-Gordon model. The strong-coupling phase is to the left 
of the KT-transition line 5*+ , the weak-coupling phase to the right of S+ . 



In spite of the strong-coupling divergence in II and III, the perturbative scaling equations allow us to analyze 



the size of the mass gap by identifying the mass M with the scaling invariant of (21). One e.g. finds the following 
expressions for small u > 



MocA for »f (22) 




for(l-^-u) >I (along S"-) (23) 



for ( I - £- + u ) > I (along S+) (24) 

In this manner one can obtain information about the renormalized low-energy scale even in the strong-coupling 
phase. But the perturbative scaling approach does by itself not lead to an understanding of the physical behavior 
associated with this low-energy scale. This situation is typical for other strong-coupling problems as well, the Kondo 
model being the paradigm in condensed matter theory llj. In combination with mappings to other exactly solvable 
models these shortcomings can sometimes be partially overcome, see Sect. II. C below. However, in general there is 
considerable interest in theoretical methods that can solve strong-coupling problems in a controlled way. Therefore 
the flow equation approach might be an interesting tool also for other strong-coupling problems by removing some of 
the above shortcomings. 

The phase diagram Fig. [l] remains essentially unchanged in higher loop orders [T^ . For latter comparison with 
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the flow equation solution it is interesting to also write down the 3-loop result for the mass gap on S- (that is for 
(3^ = 8n{l - u)) in the limit u^O 

M cx Au^/^cxp ^-^^ . (25) 

Notice the it^/^-prefactor as compared to the 2-loop result (p3|). Higher loop orders beyond 3-loop should only affect 
the proportionality factor in (^) . 



C. Integrable structure and relation to other models 

The sine-Gordon model is one of the best studied integrable models, which makes it a very suitable test model for 
our new approach. Its spectrum was obtained exactly from an inverse scattering solution JT^ and its S'-matrix was 
calculated by Zamolodchikov For a recent review see Ref. ||l9[| . 

In the strong-coupling phase the exact solution confirms the scaling behavior (|2^) for the mass M of the solitons 
and antisolitons. The exact S'-matrix jl^ also shows that for small rapidity differences these solitons and antisolitons 
behave as fermions. This important observation of a change in statistics for the low-energy excitations will be 
reproduced in our flow equation framework. 

In addition, the exact solution shows that new features appear for < 4tt: Soliton-anitsoliton bound states 
(breathers) emerge in the spectrum with excitation energies smaller than 2M. There is one breather for Stt/B < 0^ < 
47r, two breathers for 27r < /S^ < Stt/B etc. 

The sine-Gordon model is related to other integrable models like the spin- 1/2 X-Y-Z chain and the Thirring model. 
Since in particular the relation to the Thirring model will be fruitful in the sequel, it is useful to sum up some of its 
main properties here: The massive Id Thirring model is defined by the Lagrangian density 

Cty. = f V5*7p5'^ V' ~\9 j^A - Mi:^ (26) 



^(-) - ( ) (27) 



in terms of the two-component spinor ip{x) 

with components obeying fermionic anticommutation relations 

{^^.(x),4(y)} = %,5(x - y) , {^^.(x),^^.,(y)} = {V'](x),4(y)} = . (28) 

The current is defined by j'^ -ipjf^il; with -0 '== ip^jo and the 7-matrices are explicitly given by 

70= \) and 71= L? l) ■ (29) 



The exact Bcthe ansatz solution of the Thirring model was obtained by Bergknoff and Thacker pO|| . 

Coleman pl| has shown the equivalence of the sine-Gordon model with the Thirring model (p6D order by order in 
perturbation theory with the following mapping between the coupling constants 

1 

47r I + g/n 

One notices that = 47r is a special point of the sine-Gordon model since it corresponds to a noninteracting massive 
Thirring model {g = 0). For = An the elementary excitations of the sine-Gordon model are therefore fermionic 
with the dispersion relation ztEk with 



Ek = Vk^ + ■ (31) 
The explicit relation between the bosonic field of the sine-Gordon model and the fermionic field of the Thirring model 



was found by Mandelstam |22|. For = Att one has explicitly 
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i',{x) = ^V,{-l;x) (32) 



with Vj{a;x) from (|14|): Notice that the Vj(±l;x) obey anticommutation relations (A20) in the hmit a ^ 0. These 
Thirring fermion^ correspond to the quantized sohton solutions of the sine-Gordon model . 

The perturbative scaling approach does not "know" about the special point = An where the sine-Gordon 
model is trivially diagonalizable by using the equivalence to the quadratic Thirring model: The strong-coupling 
scaling trajectories in Fig. |^ go right through the line = 47r. The subsequent strong-coupling divergence of the 
running coupling constant is then due to the fact that one has generated the nonvanishing energy scale M. A similar 
scenario occurs in the Kondo model: There the diagonal Hamiltonian corresponds to the Toulouse point and the 
nonvanishing energy scale is the Kondo temperature Tk- 

A standard approach to avoid the strong-coupling divergence is to scale the model to the exactly solvable line = 
Att: One stops the scaling once /3^(Acff) — 47r and obtains the mass from the value of the running coupling constant. 
With this approach it is also plausible that the low-energy single-particle/hole excitations in the strong-coupling phase 
are fermionic with a mass set by the scaling invariant. Though very useful, this is an uncontrolled approximation 
since the running coupling constant is already large when /3^(Aoff) — 47r. It is therefore difficult/impossible to learn 
something about the crossover from weak-coupling to strong-coupling or about the effect of irrelevant operators at 
the strong-coupling fixed point. These shortcomings make it desirable to develop our new method that will allow a 
systematic expansion describing the full crossover flow. 

The sine-Gordon model is also related to a variety of other models like the 2c? Coulomb gas with temperature 
T — f3~'^ and fugacity z cx m, or a Id electron gas with backward scattering. For an overview of these and other 
relations see Ref. p6|| . As a final remark we also want to mention that the mapping to the Id electron gas gives a natural 
interpretation to the separatrices S± in Fig. |l| since they correspond to an electron gas with SU(2) spin-symmetric 
interactions |p6|. The sine-Gordon model with = 87r(l ± u), |u| ^ 1 therefore carries a hidden SU(2)-symmetry. 



III. FLOW EQUATION APPROACH 



A. General concepts 



The idea to apply a sequence of infinitesimal unitary transformations to a Hamiltonian in order to make it more 
diagonal has been independently put forward by Wegner ||] and Glazek and Wilson js). Wegner's original work 
focussed on diagonalizing many-particle Hamiltonians, whereas the focus in the work of Glazek and Wilson was to 
construct effective low-energy Hamiltonians for strong-coupling field theories: Such effective Hamiltonians can then 
be analyzed by standard techniques in order to find the bound state spectrum, which in turn could be interpreted as 
e.g. hadrons or mesons. Though the outlook of these approaches is somehow different, the concepts are very similar. 
In this paper we will follow Wegner's methodology. 

The main idea of Wegner's flow equations is to generate a one-parameter family of unitarily equivalent Hamiltonians 
H{B) labelled by a flow parameter This is achieved by solving a differential equation 

^^ = [r7(S),H(B)] (33) 

with some anti-Hermitian generator ri{B) = —rii^{B) where H{B = 0) = iJ is the initial Hamiltonian. One wants to 
choose rj{B) such that H{B) becomes more diagonal as B oo: Splitting up H{B) in its diagonal and interaction 
parts 

H{B)^Ho{B) + Hir,tiB) , (34) 

this amounts to requiring _ffi,it(S) becomes (in some sense) smaller for i? — > oo. In order to achieve this, Wegner 
proposed the following generator [g| 



'^Notice, however, that Vi(±l;a;) commutes with V2(±l;a;) instead of anticommutirig. "Proper" fermions can easily be defined 
with an additional Jordan-Wigner phase factor, but nothing new can be learned from this Jordan-Wigner construction. 

^The flow parameter has been denoted by £ in most other works on flow equations. In order to avoid confusion with the 
common notation where £ is the logarithm of the change in length scale in RG-equations, B instead of £ is used in this work. 
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(35) 



With this choice of ri{B) one can show 

-^TtHI,{B)<0 (36) 

and in this sense the operator Hint{B) becomes smaller along the flow. Notice that B has the dimension of (Energy)"^ 
with this choice. However, for a many-particle Hamiltonian Eq. ( |36| ) is usually not well-defined since the trace is 
typically infinite. Also higher and higher order interactions are successively generated by the system of equations ( |3^ ) 
and (|35|), which have to be truncated in some way making rigorous statements difficult. 

Still one finds that (^5|) is generally a suitable choice for achieving our goal to make the initial Hamiltonian diagonal if 
Hint{B = 0) can be viewed as a small perturbation term: Truncating the system of higher order interactions produced 
by ( |33| ) and (|3^) in some order of the coupling constant then amounts to generating a perturbation expansion in a 
renormalized coupling constant. From this point of view the flow equation approach is similar to perturbative RG. 
Matrix elements of Hi^t {B = 0) that couple states with large energy differences are eliminated in the initial stages of 
the flow (for small _B), and matrix elements coupling more degenerate states are eliminated later. This is reminiscent 
of the energy scale separation underlying the renormalization group approach, which is the suitable perturbation 
expansion for systems with largely varying energy scales. 

Explicit applications of these ideas have been discussed for various model Hamiltonians like the n-orbital model 
dissipative quantum systems 1^, s y stems with electron-phonon coupling ||^,^ and various other models in 
condensed matter physics ]l^,|6|,p^U9|, |2^ , p5t . In the present paper it will be shown how this method can be used for 
the sine-Gordon model as a genuine strong-coupling many-body Hamiltonian. 



B. Generator rj 



The aim of this work is to diagonalize the sine-Gordon model (^^ using the method of inflnitesimal unitary trans- 
formations outlined above. We split up the sine-Gordon Hamiltonian H{B) into a free part Hq and the interaction 
part Hint{B) 



H{B)^Ho+Hir,t{B) (37) 



with 



2N 



= J2p (ai(p)ai(-p) + a2i-p)<J2{p)) (38) 

Hint{B) = J dxdyu{B-y) {Vi{a-x)V2{-a;x~y) + \i.c.) . (39) 

In order to avoid confusion the initial parameters u and (3 in (|l^) will from now on be denoted as uq and /3o. The 
initial condition for ( ^9| ) then reads 

/3o 



Ait 

Notice that we have already allowed for a general nonlocal interaction u{B] y) in ( ^9| ) since the initially local interaction 
( ^ ) will become nonlocal along the flow (see below). Next we have to evaluate (|35|). The following commutator is 
useful 

[aj{p),Vj,{a;x)] = a exp (-^\p\ + ipx^ Vj,{a;x) (41) 

leading to 
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[J2p<^2{-pW2ip),V2{a;x)] = -id^V2ia;x) . (42) 

p>0 



Thus wc find the foUowing generator 



-2i I dxdy^^^^^{Vi{a;x)V2{-a;x-y)+h.c.) . (43) 



C. Commutator [i],Ho] 

To study the flow generated by 77 we first look at the commutator [rj, Hq]- Using (^) one easily shows 

[r,,Ho]=4. J dxdy^^^^^^{Vi{a;x)V2i-a;x-y) + h.c.) . (44) 
Comparison of the coefficients on the left-hand side of (|3|) with (E3) gives 



du{B;y) _ dMB;y) 

dB dy^ ' ^ ' 

where possible contributions from [77, Hint] are still missing. Eq. (|45| ) has the character of a diffusion equation: The 
initially local interaction becomes increasingly non-local along the flow. In terms of Fourier coefficients 

u{B;y) = J2^^iB;p)e-'P'' (46) 
p 

one finds the solution 



47r2a2 V L , 

One sees explicitly that matrix elements u(B;p) coupling states with large energy differences \p\ are eliminated in the 
early stages of the flow (for small B), whereas matrix elements coupling more degenerate states are decoupled later 
during the flow. This is a generic feature of Wegner's generator (|35|). 

Later we will see that (|47| ) is modified due to higher-order contributions. Therefore we introduce a more general 
parametrization 



v{B-p) (48) 



with a running coupling u{B), initially u{B 0) = uq. The differential equation for the coefficients v{B;p) now reads 

dv{B;p) 



dB 

with the initial condition v{B — 0',p) — 1. 



^-Ap-'viB-^p) (49) 



D. Commutator [i],Hint] 
1. General properties 

The evaluation of [77, i?int] is the key calculation in the flow equation approach. We first look at some general properties 
of such commutators. Let Ai, A2, Bi, B2 be arbitrary operators with [Aj,Bji] = 0. We define *0* as the operator 
with its ground state expectation value subtracted 
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(50) 



with the notation (O) = (i7|0|f2). One easily shows 

[AiSi, = (AiA2)(SiS2) - (A2Ai)(S2Si) (51) 

+ (i?li?2) * ^1^2 * - {B2B1) * A2A1 * 
+ (^1^2) * B1B2 * -{A2A1) * B2B1 * 

+R 

with 

R = *AiA2 * *BiB2 * — * A2A1 * *B2Bi * . (52) 

In general R leads to the generation of higher order interaction terms during the flow. R vanishes if the operators 
fulfill the following exchange relations 

A1A2 + e'*A2Ai = c 
B1B2 + e-"l'B2Bi = c (53) 

with fixed cf) and c. E.g. for = these are fermionic anticommutation relations, or for (f> = tt bosonic commutation 
relations. Then no higher-order interactions are generated and it is possible to close the flow equations without 
approximations. For general Po in the interaction term we will, however, have to develop a suitable approximation 
for R in the next section. 



2. [ri, Hint] in the sine-Gordon model 

There are two structurally different commutators of vertex operators generated by [77, i?int] in the sine-Gordon model: 
[Vi V^,V^V2 ] and [V^ k/, ^1 ^2] (or equivalently [V^V^ ,¥^¥2]). The first term 

[Vi{a;xi)V2i-a;xi - yi), Vi (-a; 2:2)^2(0;; a;2 - ^2)] (54) 

and its hermitian conjugate will turn out to be the leading contributions and are discussed first. Eqs. ( A16| ) and 
(|A17|) give 



{Vi{a;xi)Vi{-a;x2)) = s-^" , {Vi{~a;x2)Vi{a;xi)) ^ s'^ 

2 2 

{V2{-a;xi ~ yi)V2{a;x2 - y2)) = S2"' , {V2{a;x2 - y2)V2{-a;xi ~ yi)) ^ (55) 



with 



27r 

si = — (i{x2 - xi) + a) 
2tt 

S2 = — {i{xi -yi- X2+ 2/2) + a) ■ (56) 

Using ( ^l| ) we then find 

[Vi{a;xi)V2{-a;xi - yi),Vi{-a; X2)V2{a; X2 -1/2)] (57) 
22 22 

— S2 ~ S2 

2 2 

*Vi{a;xi)Vi{—a-,X2)*—S2°' * Vi{—a; X2)Vi{a- xi) * 

_ 2 2 

Si" *V2{-a;xi - yi)V2{a;x2 - y2) * -Si"' * V2{a; X2 - y2)V2{-a; xi - yi) * 

+R 



ith 



R = :¥Vi{a; xi)Vi{-a; X2) * *V2{~a;xi - yi)V2{a;x2 - ^2) * 

- *Vi{~a;x2)Vi{a;xi) * *V2(a;x2 " y2)V2{-a;xi - yi) * . (58) 
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The key approximation in our method is to use an operator product expansion (OPE) in higher-order interaction 
terms Hke R , and then to neglect contributions with larger scaling dimensions (more irrelevant terms in the RG-sense). 



From (A18) we e.g. conclude 



*Fi(a;xi)l^i(-a;x2)* = sr"' | za(x2 - xi) ^ /He-tH-'P-iai(p) + . . . ] , (59) 



where we have neglected the higher order terms in (A18). Notice that the c-number contribution has already been 



removed by subtracting the ground state expectation value. Putting everything together gives 

R^a^{x2-xi){xi^yi-X2+y2){s^"'s^''" -s^'^'s^"") ^ /Me-t(lpl+kl)-P-i-9(-i-w)^^ (p)^2(g) • (60) 

The first and second term in (|57| ) are c-numbers and describe a shift in the ground state energy. This is of no 
particular interest and we will not look into it. The various other terms generated in [ri,Hint\ are discussed in the 
next subsections. 

3. R-term 

The i?-term in (|5^) leads to the following contribution from [?/, i?int] 
-2i J dxidx2dyidy2 ^^^^^^^^ u{B\y2) 

x2a^(a;2-:Ei)(xi -yi-X2+y2)(sr"'52-"' -sr"'s2-°') ^ /Me-t(lpl+kl)-P-i-«(-i-vi)^^ (p)^2(g) 



2o? 



^ ' ^|fc|ifcai(fc)a2(-fc) (61) 



with coefficients 

tk= y"d.2iciz2(iz3e-"l'=l-*'=^^^^^^^u(B;0i + Z2)^3(^2-Z3) ((«^3 + a)""'(«(^2-Z3) + a)""' -h.c.) . (62) 
Except for an (unimportant) initial transient where B a? , the zi-integral leads to the following expression 

tk = -2Tna^-^''^"^pu{B-p)u{B;^k-p) j dxe'^^+P^"'' I{x) (63) 

with 

I{x)^ j dvy[x-y)([l + iy)-''\l-iy + ix)~''" -[l-iy)-^\l + iy-ix)—"^ . (64) 
Writing 1 + i (j/ ± f ) = r±e'">'^ with 



X 



2 



i± = arcsin^e [-|,|] (65) 



leads to 



I{x) ^2ijdy{^~ 2/2^ (r-+r_)-"' sin {a\c^. - <^+)) 



(66) 
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The flow is dominated by the term decaying most slowly with B, which corresponds to the large- a; behavior of I{x). 
Therefore we can approximate 



r± 



\z±l\ 



where z = 2y/x, and find 

I{x) = 2i 



dz (1 - zJ) |1 - sin(a^(0_ - 0+)) 



(67) 



(68) 



With the above approximation one has 0± £ §}• Thus the only contributions to I{x) come from regions with 

0- 7^<^+ 



= -7r => x>0, -l<z<l 

>,= TT => X<0, -1<Z<1 



(69) 



leading to 



I{x) = — 2isin(Q! 7r)sgn(2:) 

^3/2 

= 2i ■ 



dz{l - Z?f 



r(a2_l)r(|-a2) 



sgn(a;) 



With (|8|) and (||) the sum over p in (|3|) gives 



2q 



ax 



ik 

8B \16B ' y 



exp 



a^x^ ikax 



32B 



2Bk- 



(70) 



(71) 



The final step is to perform the a;-integration in (|63|). This can be done in closed form leading to hypergeometric 
functions. However, the flow is determined by the IR-limit fc — )■ where the integral is simpler 



r(a2 - 1)' 



For the full fc-dependence we write 



tk=tk=of{a';kVB) 



(72) 



(73) 



To leading order the only information that we will need about f{a'^\x) is that it falls of rapidly to zero for large 
arguments \x\ ^ 1. For example one easily shows 



/(q,2 ^ l-x) = e'^^\l-%x') 

f{a^ ^2;x) = e"^'-"^ - \/2^ x erf (\/2a;) e"^^' 

These functions are depicted in Fig. ||. 

Putting everything together the i?-terms in ( p7| ) from [77, iJjnt] contribute 



(74) 



a 



32 /32B\ 



\ J 2r(a2 - 1) 



-7i2(B) Y,\k\f{a^;kVB)a,{k)<j2{-k) 



(75) 



An important observation can be made for a — 1: Due to the divergent F-function in the denominator, the term (|75| ) 
vanishes for /3q = Att. This is an immediate consequence of the fact that in this case the vertex operators describe 
fermions. The interaction term of the sine-Gordon model is then simply a quadratic term in the fermions and no 
higher order interactions are generated during the flow, therefore according to ( |53| ) R = for all B. For /3q = At: 
we will be able to solve the flow equations without any approximations, thereby recovering the equivalence to the 
noninteracting Thirring model discussed in Sect. II. C. This demonstrates a fundamental difference of our approach 
to perturbative scahng, where the scaling trajectories go right through the line /J^ — An (see Fig. ft). 
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FIG. 2. f{a'^;x) for = 2 and = 1, see Eq. @. 



4- -^diag 

Let us next look at the fifth and sixth term in (|57|). The total contribution from [?7,i7int] to terms of this structure is 

(76) 



[ri,Hint] > -2i / dxidx2dyidy2 ^^'^ m(-B; 7/2) 

J dyi 



x|si" *V2{-a]Xi-yi)V2{a;x2~y2)*~s^" * ^2(0;; X2 - ?/2)V2(-a; a;i - yi) * 

2 2 

+3^" * T^2(a;a;i - ?/i)F2(-a;a;2 - 2/2) * -s^" * V'2(-q;; a;2 - y2)V'2(a; a^i - 2/1) * 

For simplicity we will only look at the first term in this expression since all the other terms can be treated likewise. 
We first exchange the two vertex operators using (A22) as this will lead to a normal-ordered expression below 



/du(B' yi) 2 3°' 

dxidx2dyidy2 Tr "(-S; 2/2)57" * V2(a; 2:2 - j/2)V"2(-Q:; Xi - 2/1)* 
dyi 

with si, 32 from (|5^). We can rewrite this in terms of Fourier transforms (a > 0) 

p 

substitute ^2^X2+ xi and perform the integral over xi. This leads to 
L 



(77) 



(78) 



27r 



dx2dyidy2 u[B; 2/2) 



dyi 



^ *V2(a; fc)V2(-a; fc) * {1x2+0) 



[i(a;2 -2/2 + 2/i) + a]" 



[-1(2:2 - 2/2 + 2/1) + a]" 

dxdyy^pu\B-p)e-'''-+'Py ^V2ia;k)V2{-a;k) ^ [z(x + ;/) + a]'"' /'"^ + , , (79) 



k,p 
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where we have employed (|4q). Next the y-integration can be done using 



dy {iy + ix + ay' e'Py = e-'P'^lpl" 



27r 



(80) 



which is vahd in the hmit \ap\ <^ 1: This holds except for an (unimportant) initial transient with wavevectors of order 
a~^. We find 



r(a2) 

Next we can approximate 



^) / f^a; X!Z!^"'"^(^'P)^~*^''^^^''*^2(a;fc)V^2(-a;fc)* 



fe p>0 



\ix + a]" 
—ix + al" 



— ?a; + al° 



cos(7ra^) + i sin(7ra^) sgn(a:) 



(81) 



(82) 



using the same reasoning as in (A24). This gives 
IGtt^ ( L 



27rcos(W)^A:" u^{B;k) * V2{a; -k)V2{~a; -k) * 

k>Q 

+isin(7ra2)^^p"\2(B;p) f dxe-'^^+P'^^sgni^x) *V2{a:-k)V2{-a;-k) 



k p>0 



2 

2«"'''S^(^)" (cos{na')J2k'^\^B;k) , V2{a; -k)V2{-a; -k) * 



+ -sm{TTa^)y y p^\^{B;p)^— * V2{a; -k)V2{-a; -k) 



(83) 



In order to do the sum over p in the second term it will be sufficient to use the approximate solution v{B\p) — e""*^^ 
from (^9|). Deviations from this approximate solution essentially only occur close to the strong-coupling fixed point 
= 1, where the second term vanishes anyway. This p-summation leads to an integral of the type 



h{a^\ x) = P / dy e 



y-x 



ie-"'|a;r'Re<!i"' 



r 1 



r ( - Y^-a^ ) - «sgn(a;)r 



I -a' 



(84) 



where T(s,z) denotes the incomplete F-function. One easily shows h{a^;x = 0) = r(Q!2/2)/2 and the asymptotic 
behavior for I a; I ^ 1 



h{a^] x) 



2x 



+ 0(x-2) 



(85) 



with a smooth crossover in between. 

It will be convenient to use the normalized operators 



dcf 



Sj{a] k) 
S]{a-k) = 



-, 1/2 



V,{-a;k) 



(86) 



l-a 



2-1 1/2 



Vjia;k) 



with the properties (see (|A7|)) 
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Sl{a;-k)\n) = S^{a;k)\n) = Sl{a;k)\n) = S2{a;~k)\n) =0 Vfc > 
and the normalization 

{Si{a;k)Sl{a;k')) = (5|(a; fc)S'2(a; fc')) = Skk'0{k) L/2tt 
{Sl{a;k)S^(a;k')) = {S^ia; k)Sl{a; k')) = Skk'0{-k) L/2tt 

for \ak\, \ak'\ <C 1 as follows easily from ( |A16| ) and ( A17 ). We express in terms of these operators and find 
X (cos{na')J2{akf"'-\'{B;k)Sl{a;-k)S,{a;-k) 

' ^ fc>0 

+ - sin(7ra2) y^{ak)'^^-^{8B/a^)'°'^^^h{a^; VSBk) Sl{a; -k)S^{a; -k) 



(87) 



/c>0 



- sin(7ra2) ^ {ak^^ -^{%B j a^y^^ l'^h{a^; -VSBk) * ^^(a; k)S2{a; k) 

TV ' ^ 



k>0 



(89) 



In the first two terms we do not need the subtraction operation * * anymore since the vacuum is alr eady annihilated 
by them. The third term does not yet annihilate the vacuum. This can be easily achieved by using (A2f:). However, 
already the second term will turn out to have hardly any effect, and the third term is again smaller than the second 
term. In order to simplify our expressions we therefore omit the third term in the sequel, although there would be no 
problem at all in carrying it along as well. Let us now also collect the other terms from ([Tq) leading to 



Y ( cos{na^) {akf^'-^v^iB; k) + - 8m{Tra^){ak)'''-\8B/aY"' ^^Ha^^VSBk) 



fc>0 



(^Si{a; —k)Sl{a; — fc) + 5*1 (a; k)Si{a; k) + 53(0:; ~k)S2{a', ~k) + 6*2(0;; k)Sl{a; k) 



(90) 



Since a generically flows as a function of _B, this implies that vertex operators with different scaling dimensions 
contribute to each wavevector k. However, most of the contribution to a given fc-vector occurs in a narrow range of 
the flow: We will see in Sect. IV. C that for a given k the main contribution occurs when B m B^ with 



Bk 
Bk 



dof 1 

dof 1 



AkVWTW 



(weak-coupling phase) 

(strong-coupling phase with mass gap M) 



(91) 



In order to simplify our notation we therefore use a single scaling dimension a corresponding to each k with a — a{Bk)^ 
The newly generated term in the Hamiltonian can then be written as 



k>0 



with 



(fc) 1^' {a{Bk); k) , P}{k) - S]{a{Bk);k) 



(92) 



(93) 



and according to (|90| ) 



doj{B;k) 
dB 



Au^jB) 
r2(a2) 



cos(7ra2) (aA:)2"'~iw2(B; k) + - sm{TTa^){ak)'''-\8B/a^)-'''/^h{a^;V8Bk)] (94) 



with ijj{B — 0;k) — 0. Using Eqs. ( p2| ) one can easily check [Ho, Hdiag{B)] = 0, therefore Hdiag{B) can be interpreted 
as diagonal. 



This approximation becomes exact in the low-energy limit, e.g. in the strong-coupling phase for |fc| <C M. 



16 



5. H, 



So far we have not discussed the commutators with the structure [V\ V2 ,Vi V2] and [V\ , ] that are also 
generated by [77, i/int]- From (^) one concludes that these commutators contain only the i?-term. The operator 
product expansion in the i?-term then generates interactions with the structure Vi(2a)V2(— 2a) etc., that is only 
terms with larger scaling dimensions. These interactions are neglected in the present approximation, just like the 
higher-order terms in (|5S|). We formally sum up these neglected terms with larger scaling dimensions in iJjcs- 



E. Commutator [7;, Hdim 



We also have to study the effect of the infinitesimal unitary transformations on the newly generated terms (^2|). An 
overlap exists essentially only for wavevectors of order B~^/^. For notational simplicity we can therefore use the 
running scaling dimension a = a{B) in (B2|) and arrive at the follwing commutator 



(95) 



dy 



k>0 



X ^[Vi(a; x)V2{~a; x — y), Si{a; —k)Sl{a; — fc) + S'|(a; k)Si{a; k) 

+S\{a\ —k)S2{a; —k) + 52(0; k)S\{a] k)] + h.i 

A typical contribution comes e.g. from 

[Vi(a; x), S\{a] k)S^{a; k)] — Vi{a; x)S\{a; k)S^{a; k) — S\{a] k) * S^{a; k)Vi{a; x) * —S\{a; k){S^{a] k)Vi{a; x)) . 

(96) 

The first and the second term on the rhs lead to normal-ordered interactions with larger scaling dimensions and are 
therefore neglected (or formally contained in iJrcs)- The third term on the rhs gives rise to interactions of the type 
-ffint leading to 

Vl, ^^diag] » -2i f dx dy ^^^^] ^ uj{B; k) 



dy 



k>Q 



X (^i{a; x)S2{a; -k){V2{-a-x -~ y)S\{a.; -k)) - Vi{a; x)S2{a; fc)(S'](a; k)V2{~a; x - y)) 
+ {Vi{a;x)Sj^{a; ~k))Sl{a; -k)V2{-a;x - y) - {S^{a;k)Vi{a;x))S\{a; k)V2{-a]x- y) 
+h.c.^ 

-2. fdxdy'-^ Y.^iB;k) ^^'^'^ 



fe>o 



L \ 27r 



Jk{x-y) 



X (vi {a; x)S2 {a; -fc)e-*'=(^-J^) - Vi {a; x)S2 {a; k)e' 

+Sl(a; -k)V2{-a; x - y)e'^'= - S'J(a; k)V2{-a- x - y)e-'^'' + h.c}j 

-An / dx^uj{B;k)ku{B;k) 



r(a^)2^ 



L V 27r 



X [Vi (a; x)S2 (a; -fc)e-''=^ + Vi (a; x)S2 (a; k)e'^'' 
+51 (a; -k)V2{-a\ xY^"" + ^{(a; k)V2{-0i\ x)e-'^'' + h.c.^ 
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-8TrLj2^iB;\k\) \k\ uiB; k)- 



1 fLlk 



T(a2) V 27r 
This term has to be compared with ( |39| ) 

-ffint = / dxdyu{B;y) {Vi{a] x)V2{-a; x - y) + h.c.) 



Sl{a;k)S2{a;k) + Sl{a;k)S^{a;k)^ . (97) 



2nLj2n{B;k)-l 



m 

r(a2) V Svr 



fc)S'2(a; fc) + S'2(a; fc)S'j^(a; fc) 



Eq. ( p7| ) therefore gives an additional contribution to the previous flow equation (^) and one finds 

dv{B; k) 



dB 



-Ak^viB; k) - A\k\ uj{B- \k\)v{B; k) 



(98) 



(99) 



We have seen that the term (|75| ) 



F. Flow of 



1. Unitary transformation e 



Y,MB)\k\a,{k)a2{~k) dB 



with 



(100) 



(101) 



is generated during the flow. This term is not contained in the original sine-Gordon Hamiltonian. We will now show 
how this term can be eliminated by an additional unitary transformation with 



U ^^tpp {ai{p)a2{~p) - cti(-p)ct2(p)) 

p>0 



(102) 



with suitable parameters ipp [p6| . 

Let us first write down some general properties of this unitary transformation for general ipp. The bosonic fields 
are transformed according to 



e ^cTi {p)e^ — (Ti (p) cosh tpp + a2 (p) sinh il)p 
e~^(T2 {p)e^ — (72 [p) cosh 7/;^ + ai [p) sinh ipp 



and vertex operators as 

e^^Vi{a;x)e" = : exp 

e~^V2{a;x)e^ = e*" : exp 
where 



a ^ cosh ^p e 2 
—a cosh V'p 



CTl(p) 



: : exp 



a ^ sinh Vp^^ e-t IPl-^P^(T2 (p) 



(103) 



(104) 



\\p\-ipx 



<72{P) 



p#0 



exp 



-a^sinhV'p^^e-tlp|-'P^CTi(p) 



C = — sinh V'p sinh i/j^p ■ 



(105) 



p>0 
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2. e-^He^ 



We take the point of view that ( 100 ) has been generated infinit esim ally by integrating the flow equations from B to 
B + dB. We therefore apply the above unitary transformation ( |102| ) to II 

(106) 



n{B + dB) — > e-^H{B + dB)e^ 
Let us first investigate the effect on Hq: For the choice 



-^dB 



(107) 



the transformation Hqc^ reproduces i?o a nd g enerates an additional term that annihilates ( |100[ ). This can be 
shown easily by using the transformation rules (103). Notice that terms of order -0^ and higher can be neglected since 
ijjp is of order dB. 

Next we have to find the effect of this transformation on Hint{B). Using (104) one finds 



e ^Vi{a;x)V2{—a;y)e^ — Vi{a;x):exp 



p#0 



xV2{-a;y) : exp 



V2(p) 



(108) 



where we have again used that tpp is of order dB. Using an OPE, we can combine the first two terms into a vertex 
operator with a modified scaling dimension, and likewise for the second two terms. The calculation proceeds along 
similar lines as in (A18): 



Vi{a; x) : exp 



a 



, p 



-ipy 



exp 



X exp 



x-v\p\ - 



p>0 



exp 



^ , vH -f 



Pl-'f^cri(p) 

"lf|-'P%l(p) 



a 



0-1 b) 



p>0 



exp 



p<0 



X exp 



X exp 



p>Q ^ 



p<0 



where 



p>0 



exp(— ap — ip(x — y)) 



(109) 



(110) 



From (74) we know that 4'p is^i^s of rapidly on an energy scale of order B ^1"^ . The leading behavior of the sum is 
therefore (except for an uninteresting initial transient) 



p>0 



exp(— ap — i'pix — y)) 



leading to (compare Eq. (A13)) 



C _ f 27r(a + iix - y)) 



a i/'p=o 



(111) 



(112) 
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Here a = s\/B takes into account the UV-cutofF generated by the decay of the functions /(a^; x) from Fig. ^. The 
actual value of the proportionality constant s will only affect our results in next to leading order as will be shown 
later. Still we give its value here for latter comparison with the three loop scaling results: For = 2 one finds in a 
somehow lengthy calculation 

ln,s = ln2+ J-^ , (113) 
where 7 ~ 0.577 is Euler's constant, s depends only weakly on a^, therefore we will use this value throughout. Next 



exp 



a 



E 

p>0 



p 



e 2 



bl-ip^(l + ^pe-*P(^-^))CTi(p) 



exp 



q;(1 + Vp=o)X! 



e 2 



\p\—ipx 



0-1 (p) 



p>Q 



X exp 



iix - y)aijp^o E /H e-tif l-PVi(p) + 0((x - yf) 

p>0 



^Viia{l + ,Pp^o);x) X l + z(x-y)aV'p=oE^^'^''''"'''"'^i(^')+0((^~2^)') 
— Vi(a(l + 'ipp^o);x) + more irrelevant terms 



(114) 



up to less singular (more irrelevant) terms. We formally include these more irrelevant terms (they have the structure 
of products of vertex operators multiplied by derivatives of the bosonic field) into i?rcs and neglect them from now 
on. It should also be noted that we have approximated the flow of the scaling dimension in the vertex operator 
by restricting ourselves to the flow in the IR-limit a — *■ a(l + ipp^o)- At first sight this seems to be a problematic 
approximation since the functions f{a^;x) in (fz^ ) are nontrivial. However, the effect of the unitary transformations 

is only cumulative on the low-energy scale where f{a^;x) 1. Hence it is possible to restrict ourselves to the 
IR-hmit. 

Putting everything together we find 



Hi,,t{B + dB) 



2'KS\/B \ 



(115) 



dx dy u{B; y) 



VAa{l-^^^^^dB),x\V, 



-a{l-^^^^4^dB;x-y 



h.c. 



Similar to the OPE above we have set cc — y = in (112). that is we have neglected less singular terms. 

Finally, we have to investigate the effect of the additional unitary transformation (106) on -ffdiag from 
similar to the action on the interaction term, except that here we find terms of the structure 



I). This is 



-""Vi (a; x)Vi {-a; 2/)e^ = Vi (a; x)Vi {-a; y) 



: exp 



a 



Z — ^ ^ 7) 



p#0 



X : exp 



p#0 



-ay _ipp 



CT2(P) 



(116) 



and likewise for j = 2. The third and fourth term can be combined using an OPE (A19) and one easily notices that 



the only surviving term is a constant 1 since all other terms are of order — 0{dB'^). Hence 



i/diag(S + dB) e-^i7diag(S + dB)e^ = -ffdiag(S + dB) 



(117) 

Summing up, we have looked at another infinitesimal unitary transformation ( |106| ) that acts on the Hamiltonian 
during the flow equation procedure in addition to the generator (p3|). An alternative viewpoint is to say that the full 
generator of the flow now takes the structure 
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VnoAB) = -2i Jdxdy ^^^^^ {Vi{a;x)V2{~a;x^ y) + h.c.) 

+ \Y,Wp{B) {a,{p)a2{-p) ~ a^[-p)a2{p)) . (118) 



p 



3. Flow equation for (3 



From (|TT5|) we can read of that the additional infinitesimal unitary transformation generates a flow in the scaling 
dimension of the vertex operators 



a{B) — » a{B) 



a(B) (l^^^^I^dB 



dB 



— = -Wp=o{B)a^B) 



(119) 



32 /32B\ 



a\B) 



2T{a^{B) - 1) 



u\B) . 



We can already see that = 1 is a fixed point of the flow equation approach due to the diverging F-function in the 
denominator of 
According to ( 



11£ ) . This will he one of the key results of our new approach. 

115) this flow of the scaling dimension now induces a flow of the coupling constant u{B;y) 



u{B;y) — > u{B;y) 



'27rs\/S^ 



-a^(B)u)p=o(-B) dB 



u{B; y) 1 - {B)wp^a [B] dB In 



'2■ns^/B^ 



du{B-y) ,da\ j 27rs/B 



(120) 



The solution is straightforward 



u{B; y) = u{B = 0; y) exp I / dB' ^ \n ' ^""^^ 



dB' 



.(S = 0;y)exp^ —In^—^-^ 
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= u{B^Q-y) 
Using the parametrization ([48| 



27ra\ 



Q^(B)-Q^(0) 



q2(_B)-q2(0) 



exp 



dB' n 

dB' \ a2 



(121) 



viB;p) 



(122) 



we s ee that this can be most conveniently expressed as a flow equation for the running coupling constant u{B) in 

(Hi) 



u{B) = uo 



a^(B)-a^(a) 



exp i 



dB' V 



(123) 



Introducing the dimcnsionless logarithmic flow parameter i 
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dcfl, (i2B\ 



one can show by partial integration 



1 



da" 



fi2B' 



2 ./o dB' V a 



di' a'i£') + a'{e)£ 



Using this we can sum up the resuhs of this section in the foUowing two equations: 



u{£) = uq 
da^ 



32/ 

2x Q"(i!)-a"(0) 



exp 



d£' a'^{£') + a^{l)£ 
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(124) 



(125) 



(126) 



(127) 



These two equations constitute the key results of this work. Eq. (|127| ) describes the flow of the scaling dimension 
under the flow equation procedure, and from Eq. (126) it follows how this induces the flow of the running coupling 
constant u{£). Therefore these equations will serve as a generalization of the scaling equations derived in perturbative 
renormalization theory in Sect. II. B. The value of the constant s in these equations will turn out to affect our results 
only in next to leading order. 



IV. SOLUTION OF THE FLOW EQUATIONS 

A. Summary of the flow equations 

In this section we will sum up the resu lts for the flow of the sine-Gordon Hamiltonian under the effect of the 
infinitesimal unitary transformation (118) as derived above. For general B the sequence H[B) of unitarily equivalent 
Hamiltonians takes the form 



Here 



dx 



H{B) = Ho + Hint{B) + i?diag(B) + H,,,{B) 



2 ^ ' 2\dx 



HintiB) = J dxdyu{B;y) {Vi{a{B);x)V2{~a{B)-x~y) + h.c.) 



(128) 

(129) 

(130) 
(131) 



fc>0 



with Pj{k) given by (M 



Pjik) = 



Tia^iBk)) (L\k 



2t:L 



27r 



l-a^(Sfc) 



1/2 



dx ( 



'Vj{-aiBky,x) 



(132) 



and Bk from (|9^). The operators P-{k),Pj{k) will turn out to be the soliton creation and annihilation operators. 
They are normalized according to (pq). Notice Hdiag{B)\fl) = 0. 

iJios contains the neglected terms and will from now on be omitted in our analysis. At any given B-scale these 
neglected terms have a larger scaling dimension than the interaction term HintiB): They are more irrelevant by at 
least two spatial derivatives. Notice that i?rcs vanishes for = An since then no approximations are made. 

Summing up the differential flow equations for the parameters in H{B) we have 



dv{B;k) 
dB 
d{auj{B;k)) 
d{B/a^) 



-4fc^u(S; k) - 4fcw(B; k)v{B; k) 



(133) 



Au'^jB) 
r2(a2) 



{cos{-Ka^) {akf'^^-\^{B- fc) + ^ s\n{TTa^){ak)'^^-^{^B/a^)-"'^ l'^h(a^; VSBk)^ (134) 
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for A; > 0, and v{B; —k) = v{B; k) symmetric in k. For notational convenience we have written a{B) without its 
argument B in these equations. h{a^;x) has been defined in (84) and 



u(B) (2-Ka 



V 



v{B;p)e' 



(135) 



The initial conditions are i?rcs(S = 0) = 0, uj{B = 0; fc) = and v{B — 0:k) ~ 1, whereupon (128) takes the form ( |l8| ) 
of our original sine-Gordon Hamiltonian. For ag < 1 one should in fact be more cautious and start the integration 
only at B — a?: The above equations hold only for B }^ since we have used \ak\ ^ 1 in our calculation. One can 
easily verify that there is no flow of the parameters for i? <C a^, and the simplest way to take this into account is to 
pose the initial conditions at B = a? . 

As we will see later the flo w of the parameters is such that v{B — oo; fc) — for /3q > 27r. In this parameter region 
the Hamiltonian H{B) from (128) therefore becomes diagonal with Hint{B — oo) = in the limit i? ^ oo as expected 
under the flow equation procedure 



H{B = <x) = Ho + ffdiag(S = oo) 



(136) 



notice [f/"o, -ffdiag(S 



0. 



The equa tion s ( |133| ) and (134) have to supplemented with the differential equat ion governing the flow of the scaling 
dimension ( |127| ) and the thereby induced flow of the running coupling constant ( |126| ). It is possible to rewrite (127) 
in a more conventional form for comparison with the perturbative scaling approach. Introducing a new function 



u{i) =^ Mo 



exp 



2£- 



dt' a' {I') 



(137) 



we can rewrite ( 127 ) as a set of two coupled differential equations 

a4(£) 



d£ 

de ^ 



r(a2(£) - 1) 



a' 



•■{£)) u{£) 



4-2-2^"^ 



(138) 
(139) 



with the initial conditions u{£ = 0) = u q an d a{£ = 0) = /So/y/An. For convenience we have use d the dimensionless 
flow parameter £ = ln{32B /a^) from (124) in these equations. Notice that the second term in (13£) is of order 
and does therefore not contribute to the leading behavior for small uq. Since our present flow equation expansion has 
not taken all the terms in order into account anyway, we can omit this term and arrive at 



da^ 
~d£ 
du 



aH£) 



r(«2(^)_i) 



u\£) 



-^{2^a^{£))u(£) 



The running coupling constant u{£) from ( |126|) can also be expressed as 

u{£) = u{£) exp((a2(£)-2)£) . 
In terms of the sine-Gordon parameter (3 these equations take the equivalent form 

d/3-2(£) 1 



(140) 



(141) 



d£ 47rr(- 1 + /32(^)/47r; 
P^{£) 



u\£) 



and 



d£ 



u{£) = u{£) exp 



47r 



P\£) 



i{£) 



2 £ 



(142) 



(143) 
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Finally it is of some interest to express Hint{B) directly in terms of the bosonic field (j){x) and its dual <d{x) (see 
Eq. (|ll|)). After a short calculation one finds 

i?int(B) = ^ J dxdyv{B;y) COS (^(3{B) J dec{e) ^{cl){x + e) + cj^ix - y + e) + Six + e) - e{x ^ y + e))^ . (144) 

Since v{B; y) becomes more and more nonlocal during the flow, one sees that the interaction term of the sine-Gordon 
model evolves from the original cos(/30(a;))-structure to a nonlocal interaction term with the structure 

cos {W/2){c^{x) + c^{x - y) + Q[x) - Q{x - y))) . (145) 

It is also possible to express i?diag(S) in terms of 4>{x) and 8(x), however, this expression does not lead to new 
insights. 

B. Strong— coupling phase 

1. Fixed points and phase structure 



We will now work on the explicit solution of the flow equations (133), (134) and ( 142 ). First we focus on t he sol ution 
of ( |142[ ) , since kn owled ge of the flow of j3{B) and u{B) is necessary for solving the system of equations (133) and ( |134| ) 
later on. From (lli^ one concludes that there are two possible kinds of asymptotic behavior: Either/?^ (oo) = An 
or /3^(cxd) > Stt. (t\oo) — An will turn out to be the attractive strong-coupling fixed point and values /3^(cxd) > Stt 
correspond to the gapless weak-coupling phase. These flows are depicted in Fig. ^ Notice that the fundamental 
difference from the perturbative scaling equations ( pT| ) is the F-function in the denominator of our flow equation 
for (3~'^(l). Therefore 0^ = 47r is a fixed point in our approach which will be the main difference as compared to 
perturbative RG. This does not come as a surprise since the flow of 0^ followed from higher order terms in the 
commutator [r]{B),Hii-^^{B)] in Sect. III.D. However, for 0^ — An our interaction term is quadratic if considered as 
an interaction term for Thirring fermions (see Sect. 11. C), and naturally no higher order terms can be generated in 
[r](B), Hint{B)]. Another way of saying this is that the flow in is due to approximations in the flow equation scheme 
when higher order terms are generated. No flow of can occur if the flow equations close exactly. 

In both phases ■u{i) remains finite for £ oo. On the other hand, one easily checks that u{£) diverges in the 
strong-coupling phase an d va nishes asymptotically in the weak-coupling phase. Since u{£) has so far only been 
introduced for rewriting ( |127| ), its divergence in the strong-coupling phase need not worry us here. The question 
of the true expansion parameter of our approach will be discussed below in Sect. IV. D, and we will see that this 
expansion parameter is not u{£). 



In order to analyze the phase boundaries we can expand the F-function in (142) around 0^ = Stt. In leading order 
this reproduces the perturbative scaling equations 



d£ An 
d£ V An 



du f 0\ , , 

' 2---] u . (146) 



This approximation eventually breaks down in t he s trong-coupling phase as 0^{£) flows to An: Then (|146| ) is not a 
good approximation for the true flow equation ( |142| ) anymore. Notice the sign difference from (^l|) because £ from 



(124) corresponds to — In A, therefore £ is integrated from to 00. Our ffow equation approach therefore reproduces 
the conventional two-loop scaling equations if we expand around 0^ = Stt. In this way we also reproduce the hidden 
SU(2)-symmetry of the sine-Gordon model for 0^ = 87r(l ± uq) mentioned in Sect. II. C, although our approximation 
scheme does not manifestly respect this symmetry. Besides showing the consistency of our new approach with the 
conventional perturbative RG scheme at 0^ « Stt, we also see immediately that our flow equations reproduce the 
Kosterlitz-Thouless phase diagram Fig. |^ of the sine-Gordon model established with RG: In the limit of small initial uq 



^Implicitly an assumption about the behavior of v{B; k) has been made when deriving (|l42|) in Sect. III.D. One can check 
that this assumption is self-consistently justified by analyzing the whole system of equations. 
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FIG. 3. Flow of I3'^{i) under the flow equation procedure for uq — 0.1. Observe the strong-coupling fixed point 
f3'^{£ — oo) = 4tv and the weak-coupling fixed points with I3'^{£ — oo) > 8tv. 



and /3g > 87r(l + uq + 0(uq)) we flow to a weak-coupling fixed point, for /3q < 87r(l + uq + 0(m§)) to the strong- 
coupling fixed point /3^(oo) = 4tt. These two phases are again separated by a Kosterlitz-Thouless type transition 
along Pq = 87r(l + uo + 0{uq)). For the rest of this section we will focus on the strong-coupling phase. We will return 
to the weak-coupling phase in Sect. IV. C. 



One advantage of the flow equation scheme is that we can easily analyze the behavior of our model with the final 
diagonal Hamiltonian H{B = oo). However, a simpler kind of analysis is possible by identifying a low-energy effective 
Hamiltonian and analyzing this effective Hamiltonian. This will be done in this subsection. Our results will be 
confirmed by the analysis of H{B ~ oo) in Sect. IV. C later on, but the identification of a low-energy effective 
Hamiltonian allows us to make contact with the conventional scaling picture, which is very useful for providing a 
simple coherent description of the flow equation approach in the strong-coupling phase. 

Let us look at H{B) for large B such that \0^{B) — 47r| <C 1. Then we can approximately set a{B) ~ 1 in Hint{B) 



2. Low-energy effective Hamiltonian 



from (|T3^) and rewrite ( p^ 



H{B)=H,ff{B)+H^i,g{B) , 



(147) 



where 



Hcs 




Now for a = 



1 the creation and annihilation operators Pj{k), Pj (k) from ( p3| ) obey fermionic anticommutation 



relations ( |A20|) 



{P^{k),P,ik')}^S,,,^ , {P^ik),P^ik')} = {P,{k),P,{k')} = 0. (149) 



25 



One could also rewrite the kinetic term Hq in terms of these fermions and would then arrive at a noninteracting 
Thirring model ( p6|) with a nonlocal mass term as the low-energy effective Hamiltonian in the strong-coupling phase. 
However, we can also analyze the spectrum of our low-energy effective Hamiltonian Hes{B) directly by working out 
the following commutators 



[HMB),Pl{k)]^kPl{k) 
[HMB),PUk)]^~kPUk) 
leading to the dispersion relation 



u{B) 



v{B-k)Pl{k) 



u{B) 



v{B-k)Pl{k) 



Ek = \ k^ 



v{B;k) 



(150) 



(151) 



This dispersion relation describes the single-particle/hole excitation spectrum of the full Hamiltonian H{B) for 
momenta |fc| <C 1/^/B: According to (134) the terms in H^i^^{B) corresponding to such momenta are only generated 
for even larger B, therefore we can neglect the effect of Hdiasi B) f or excitations with |fc| ^ 1/Vb. In this limit we 
also find the initial value v{B; fc) = 1 unchanged according to (|133| ). Summing up, in the limit fc — » the dispersion 
relation for single-particle/hole excitations in the strong-coupling phase has the form ±Ek with 



Ek = \/fc2 + M2 



and the mass 



M 



u{£ — oo) 



(152) 



(153) 



Eq. (152) is also the form expected from exact methods using integrability |^. Equation ( 153 ) is a key result in this 
work since it describes the relation between the running coupling and the generated mass term in the strong-coupling 
regime. 

We observe that the finiteness of the running coupling u(£) in (148) in the limit ^ — > oo is of fun damental importance 
in the flow equation scheme since u{oo) sets the generated mass gap in the spectrum. Via (|143|) we can also establish 



the following expression for the "usual" running coupling u{£) in the language of the scaling equations (142). In the 
limit of ^ — !■ oo one finds 



u{£) = aMe^ 



(154) 



Since Ab oc B plays the role of an effective UV-cutoff generated by the fiow equations, this means that the 
dimensionless parameter u{Ab) diverges simply as the mass M divided by the effective cutoff Ab 



, , M 
u{B) cx — 



(155) 



which allows a simple physical picture of the diverging u{B) in the strong-coupling phase. Again, this does not imply 
the breakdown of our approach since u{B) is not the expansion parameter in our approach (see Sect. IV. D below). 



3. Scaling behavior of the mass gap 



We will now analyze the behavior of the single-particle/hole mass M in various regimes and compare this with results 
obtained with other metho ds. F irst we concentrate on the scaling behavior in the limit uo as this can be derived 
analytically: According to (154) it amounts to finding the scaling invariant 



I{£)^e-'f(u{£),p\£)) 



with some suitable function /(u,/?^) such that 



dl{£) 
d£ 



= 



(156) 



(157) 
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along the flow generated by (142). Like in the conventional scaling analysis the mass M is then given by 



Mcx/K,/3o)/« • 

One finds the same behavior as in two-loop order in Sect. II. B: For example for fixed 0^ < 8tt and uq 
easily check that = Z*'^) gives a scaling invariant up to terms in second order 

dl{i) ^ff u{e) Y\ 



1(1) X O 



(158) 
one can 

(159) 



and therefore 



(160) 



in agreement with (|2^). Likewise one also finds the scaling behavior (|2^) and (gj) since we could reproduce the 
perturbative RG-equations (146) in the vicinity of = Stt within our flow equation approach. 

Since we find agreement in two-loop order, it is of some interest to also compare with higher loop calculations 
For simplicity we focus on the ma ss ga p alon g S in Fig. |l|, that is for /3g — 87r(l — uq) in the limit uq 0. We write 
— 87r(l — v{i)) and expand ( |l38| ) and (13£) up to third order in u and v 



^ = 21*2- 4(1 + 7)m\ 



du 
di 



2uv - 4 



7 
2 



■Ini 



(161) 
(162) 



where 7 ; 
finds 



with 



0.577 is Euler's constant. It is straightforward to derive the scaling invariant from these equations and one 



(163) 



M (X Uq exp ( — - — ) /a 



l-ln8 + 7r/2 
3 



0.16 



(164) 



This should be compared with the three loop result (25) |lj] with the correct exponent t^q = 1/2. [| We see that the 
present order of our flow equation expansion is correct up to two loop order and deviates if compared with three loop 
RG. This is not surprising since we have not systematically taken all the terms in order u^{£) into account in the 
present order of our flow equation scheme and there are contributions in order m issing on the rhs of ( |162| ). 

In general no c losed analytical solution for the set of differential equations (142) could be found. Some numerical 
solutions of (142) are depicted in Fig. ^. Of course the scaling behavior (160) is reproduced by these numerical 
solutions as can be seen in Fig. ^. Finally for /3q = An one can easily prove M = ug/a exactly using the above flow 
equations as should be expected since our scheme becomes exact in this case. 



C. Properties of H(B = 00) 

Since the flow equation procedure diagonalizes the sine-Gordon Hamiltonian, we can not only learn something about 
the mass gap in the spectrum, but analyze the entire dispersion relation throughout the crossover region. The flnal 
Hamiltonian takes the structure 



H(B ^^)^Ho + i/diag(B = «3) 



(165) 



with 



The deviation from the result for r in Ref. |13| occurs because the term in order u in (L39) was neglected there 
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-Inu^ 

FIG. 4. Soliton mass M as a function of the coupling constant for various values of /3q: The full lines are constrained fits of 
the power law behavior aM oc u^'^ /^o/^t) fj.Qjjj ^jjg exact solution Ref. to the flow equation results (open circles) with the 
proportionality constant being fitted. The dashed line is the case /3o = 47r where the flow equation approach agrees trivially 
(see text). 



irdiag(B - oo) - = oo- k) (P^ {^k)Pl{-k) + Pl{k)P^ (fc) + Pl{-k)P^ i-k) + P^ {k)Pl{k) 



k>0 



(166) 
(167) 



and Pj{k) defined in (^). Notice that [iJo, i?diag(S = oo]l = and H{B = oo)\n) = 0. Using (|A26D and (|87 
(see also the reasoning below (j^)) and the normalization (js^) one finds the following single-particle/hole excitation 
spectrum of H{B ~ oo): 

• Sohton (particle) excitations with excitation energy Ef^: 

Plik)\n) for fc > 
Pj(fc)|rj) for fc < 

• Antisoliton (hole) excitations with excitation energy —Ek'. 

Pi(fc)|rj) for fc < 
P2ik)\n) for fc > 
The dispersion relation Ek is given by 

Ek^\k\+uj{B ^oo;\k\) . (168) 



In order to find uj{B = oo;A;) we next have to solve the system of equations (133) and ( 134[). A closed analytical 
solution has not been possible except for the trivial case /3q = Att where one reproduces (|3l|) exactly. However, 
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we will see below that the flow of u>{B\k) from its initial value to uj{B = co;k) occurs on the _B-scale Bk ( |9l| ) 
and is negligible for B <^ Bk or B ^ Bk- We can therefore to a good approximation replace a{B) and u{B) in 
( |134| ) by its values for B = B^ and consider them as constants. Also one can v erify numerically that the term in 
the differential equation proportional to sin(7ra^) changes the dispersion relation ( 168 ) only in relative order 1% for 
> 47r (for < 5 it e.g. affects u!{B = oo; k) less then 2%). In order to gain some first analytical insight we can 

neglect it. Notice that this approximation becomes exact in the low-energy limit in the strong-coupling phase since 



a^B) 



_B-»oo 



1. We arrive at (fc > 0) 



dv{B;k) 

dB 
duj{B;k) 

dB 



= -Ak^v{B- k) - AkLo{B; k)v{B; k) 
= 4:kv'^{B;k)ck 



with 



Cfe 



cos(7ra2(Bfe)) ii^iBk) 



\ak\ 



T^a^Bk)) a2 

Now this approximated system of flow equations can be solved easily. One finds for > > — fc^ 



j{B; k) ^ -k+ \/k^ +Ck coth 4k\/P + Ck B + arccoth 



Ck 



and for Ck > 



j{B; k) = — fc + ^/fc^ + Cfc tanh ( 4fc-\/ fc^ + Cfe B + arctanh 



In both cases one obtains 



uj{B = oo; fc) 



Ck 



(169) 
(170) 

(171) 

(172) 

(173) 
(174) 



In the case Ck < — fc^ the solution for uj{B; fc) diverges as i? — > oo. For small initial couplings uq this scenario can 
according to (171) only occu r for a^{B — ) < 1/2. This just defines our permissible range of paramters /3q > 2tt 



as mentioned above. From (172) and (173) we can also read of the justification for our above approximation in 
the system of differential equations: Nearly all the flow from llj{B = 0; fc) = to uj{B = oo; fc) occurs on the scale 
B « [4fc\/fc2~+£^]^i fa Bk with Bk from (91). One can verify numerically that for Pq > An and |fc| < 3'u(oo)/a 
the solutio n (174 ) of the approximated flow equations agrees to within 20% with the full numerical solution of the 
equations (133) and (134) and becomes exact for |fc| ^ M(oo)/a. 
Putting everything together the dispersion relation is 



Ek = \ k^ ^ cos{TTa^Bk)) 



u{Bk) 



\ak\<^^iB,y 



(175) 



.r(a2(i?,)) c 

In the low-energy limit the dispersion relation ( |175[ ) takes different forms in the weak- and strong-coupling phases: 
• Weak-coupling phase: Here a^{B = oo) > 2 and we find for |fc| ^ u(oo)/a 



Ek = |fc| 



cos(7ra^(oo)) 



1 



r(a2(oo)) 



u{oo) |afc|"'(°°)-2 



(176) 



that is a gapless spectrum with _Efe = |fc| for fc ^ 0. 
• Strong-coupling phase: Here a^{B = oo) = 1 and we find for |fc| ^ u{oo)/a 



Ek 



'fc2 



u(oo) 



(177) 



This agrees with the result (152) obtained from the effective Hamiltonian analysis in Sect. IV. B. 2, that is we 
find a gapped spectrum with the mass M = u{oo)/a. 



29 




1 2 3 
k/M 



FIG. 5. Universal scalinglimit fo r th e dispersion relation Ek and the scaling dimension of the soliton excitations obtained 



by numerical solution of (P_33) and (134): The full lines depict Ek, the dotted lines vfc^ + Af^ for comparsion, and the dashed 
lines the scaling dimension a [Bk) corresponding to the resp. wavevector. 



One can verify numerically that in the strong-coupling phase for 0^ > Air the full dispersion relation obtained by 
solving (133) and (134) is very accurately described by \/k^ + M"^ even in the crossover region: In the small coupling 
limit |tto| <C 1 there are /3o-dependent universal corrections in the crossover region k = 0{M) that vanish for /3g in 
and reach at most 3% (for f3^ — Stt).^ The respective scaling forms of the dispersion relation are depicted in Fig. ^. 

Notice that according to ( |93[ ) the scaling dimension of our single-particle/holc excitations varies continuously along 
these dispersion curves, see Figs. ^ and ^: In the strong-coupling phase one finds the initial scaling dimension a(0) 
for excitations with large momenta 3> M , and the low-energy effective Thirring fermions with a(oo) = 1 for small 
momenta |A;| <€. M. This is consistent with the exact ^-matrix results by Zamolodchikov discussed in Sect. II. C. 



/2 XYZ-chain in Ref. 
is difficult since there are 

nontrivial renormalization subtleties in mapping the continuum sine-Gordon model to the disrete spin chain (in this context 
see e.g. Ref. isl). 



Such small corrections might be expected from the exact results for the one-dimensional spin-1/ 
where the form (177) holds exactly also in the crossover region. However, a strict comparison with |27 
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FIG. 6. Universal curves for the effective scaling dimension a^{Bk) corresponding to a soliton excitation with wavevector k 
in the strong-couphng phase. The curves are from top to bottom for /3o — Sir, 7tv, Gtt, Stt, iir, Svr. 



Also notice that our elementary excitations in this section are described with respect to a transformed basis since 
H(B = oo) and H{B = 0) are related by a complicated unitary transformation. 



Finally let us look at the solution of the equations (132) and (134) for initial parameters Pq < An. We have 



seen above that the differential equation for uj{B\k) leads to divergences for fc — > if < 1/2.[| In other words 
the cos(/3o<?!'(a;))-term becomes too relevant for our flow equation approach in its present form for /3q < 27r. This is 
the reason why the parameter space of the sine-Gordon model that we can deal with in this paper is restricted to 
Pq > 27r. In the interval 2tt < /3q < An our approximations become exact in the limit /3q Att and, according to the 
observations above, eventually break down for /3q < 27r. We will investigate the accuracy of our approximations in 
more detail in the following section. 



D. Approximations and expansion parameter 



Both during the flow and for the analysis of the final diagonal Hamiltonian we have neglected the terms in Hj-cs{B) in 
(128). For any given i?-scale these terms are more irrelevant by at least two spatial derivatives than the interaction 
term iJint(-B). Our approximations are therefore similar to the expansions in renormalization approaches. To judge 
the accuracy of our approximations in more detail it is important to study the prcfactors of these terms, that is to 
find the expansion parameter of our approach. Since Hjntj B) is treated as the perturbing term, we can do this easily 
by comparing it with Hq: E.g. from the dispersion relation ( |175| ) one sees that on the momentum scale k the effect of 
Hiut{B) as compared to Hq is 



In fact the sin(7ra )-term in (134) already leads to IR-problems for Uq < 1. 



The source of this problem is related to 
the breakdown at = 1/2 and can be resolved in a likewise manner as will be shown in a subsequent publication: The 
sin(7rQ^)-terms turn out to be generally unimportant and we can safely neglect them in our present discussion also for < 1. 
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FIG. 7. Universal curves for the expansion parameter g{t) of the flow equation approach (17£). The curves are from top to 
bottom for /Jg = Stt, Ttt, Btt, Stt, 47r, 77r/2, Stt, 57r/2. Notice that (;(^) = for /Jg = 47r since then our approach is exact. 



|afc|"'(-S'=)u(Sfe)/a 



oc u{Bk) 
= u{Bk) , 



(178) 



where we have used (141). Not surprisingly the expansion parameter seems t o be the running coupling constant 
u{B) like in perturbative RG. According to the systems of differential equations ( |l42|) this is good news in the weak- 
coupling phase since u{B) decay s to zero on small energy scales {B — > oo). However, in the strong-coupling phase 
u{B) diverges according to (154). 

Still our method is a systematic approximation even in the strong-coupling phase since our solution becomes exact 
for 0^ = 47r. As we have seen in Sect. Ill, our system of flow equations closes for = 47r (a^ = 1) and no higher 
order interactions are generated during the flow. Therefore Hres{B) vanishes identically on this line. This remarkable 
observation is our main difference from perturbative RG. It is in fact even a trivial observation since according to 
Sect. II. C the sine-Gordon Hamiltonian becomes equivalent to a noninteracting Thirring model for — Att. And a 
quadratic Hamiltonian can easily be solved exactly with our scheme of unitary transformations. 

Therefore the true expansion parameter of our method is necessarily some product of u{B) and {a^{B) — 1). One 
can check explicitly that the leading terms in H^csiB) contribute like 



giB)'':^'u'iB) (a^(B)-l) 



(179) 



In this context also compare our analysis in Sect. III.D.3 where we have seen in ( |75| ) that the _R-term (that we had 
to treat approximately) vanishes linearly like (a^ — 1) as ^ 1. 

This dimensionless combination g{B) is therefore the expansion parameter of our scheme^and has to remain small 
throughout the entire flow in order to have a systematic and controllable approximation. 



^Also the combination g(B) u{B) appears, but this does not make any difference for our analysis. 
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In order to verify this let us first investigate how (q!^(_B) — 1) vanishes for _B — > oo in the strong-coupling phase: 



We define e{£) =^ a'^i^) ~ 1 and approximate ( 138 ) for large i yielding 



de{£) _ 
d£ ~ 

which can be solved easily leading to 

e{£) = e{£o) exp 



^ei£)a'M'e'' , 



[aMf 



2t _ ^2to-, 



This very fast decay for ^ — * oo avoids a divergence in our expansion parameter g{£) 



g{£) (X {aMf exp 



2£- 



{aMf 



,21 



£-.00 



(180) 



(181) 



(182) 



For the full flow of g{£) one finds universal curves in the small coupling limit wq ^ that depend on /3q. Numerical 
solutions for these universal curves are shown in Fig. |^. The fact that one finds universal nonvanishing curves for 
7^ At: means that there are nonzero corrections to our present approximations even in the limit uq 0. This is not 
surprising due to the strong-coupling nature of our model. A precise statement about the actual size of the errors is 
difficult since this depends on unknown prefactors with which g{£) enters. However, it is encouraging to observe that 
\g{£)\ remains relatively small throughout the entire flow for all /3q > An. Higher order corrections can therefore be 
expected to be small and are systematically obtainable in an expansion that takes more and more irrelevant terms 
into account in our flow equation procedure. 

Finally, for /3q < Att the maximum in \g{£)\ grows rapidly as /3q becomes smaller as can be seen in Fig. |. The 
error of our approximation therefore becomes larger as /3q 0, which agrees with the observation in Sect. IV. C that 
our present scheme actually breaks down for /3q < 2tt. This is not unexpected since the perturbing cos(/3(/)(x))-term 
becomes more and more relevant for small Pq. 



V. CONCLUSIONS 



In this paper we have developed a new approach for solving the quantum sine-Gordon model (20) 



H(i,.0)^/..(ln=w + l(|)%^c„.[,/ 



de c(e) (j){x 



(183) 



by means of infinitesimal unitary transformations as introduced by Wegner |^ and Glazek and Wilson Within 
an approximation that neglected operators with larger scaling dimensions (more irrelevant terms) we obtained a flow 
that unitarily linked the initial Hamiltonian (|183[) to a diagonal Hamiltonian 



H{B 



fe>0 



uj{B 



; k) (-fc)P/(-fc) + Pl{k)P, (fc) + PU-k)P, i-k) + P, {k)Pl{k)) . (184) 



Here the Pj{k) are soliton and antisoliton creation and annihilation operators as defined in (^3|). Their dispersion 
relation ztEk was calculated in (175). In the small coupling limit \u\ ^ 1 we found Ek = vk^ + AP (with very small 
deviations from this form) in the strong-coupling phase, and a gapless spectrum Ek = |fc| in the weak-coupling phase. 
In the strong-coupling phase, our low-energy solitons and antisolitons are fermionic (compare Fig. |^) as known from 
the exact S'-matrix solution |l3|. Within our approach their mass M can be obtained by the solution of the flow 



equations (142) 



d(3-^(£) 



1 



d£ 47rr(- l + /32(^)/47r) 
du 13^ (£y 



\£) 



d£ 



Att 



u{£) 



(185) 
(186) 
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via the relation M = e~^u{t)/a in the hmit £ ^ oo. The above equations have to be integrated from their initial 
values for £ = to = oo. Our results for the scaling behavior of the mass agree with exact methods and the 
two-loop scaling analysis (compare Fig. |^). We also reproduce the phase diagram Fig. I of the sine-Gordon model in 
our approach. 



Our equations (|185| ) and (186) are similar to the two-loop RG equations ( pi| ) except for the F-function in the 
denominator of ( |l85D that makes — 4tt an attractive strong-coupling fixed point of the flow equation method (see 
Fig. ^). This is the main difference between our approach and perturbative RG. Since the sine-Gordon model for 
= 47r can be interpreted as a noninteracting Thirring model, our flow equation procedure becomes exact at this 
point and diagonalizes the ensuing quadratic Hamiltonian easily. The expansion parameter of our approach is therefore 
not u{i) (notice that u{£) diverges in the strong-coupling phase as in perturbative scaling), but according to ( |179[ ) the 
product g{i) = u'^{i) (—1 + /3^(^)/47r). g{i) remains small throughout the entire flow for > An (compare Fig. |^). 
This allows a systematic improvement of our present approximations by successively taking terms with larger scaling 
dimensions into account in our flow equation procedure. Furthermore, it allows us to conclude that our present 
approximation already provides a good description of the crossover region, which is notoriously difficult to study 
with other techniques. For example, we worked out the dispersion relation of the single-particle/hole excitations 
for all momenta in Fig. ^. Notice that higher order terms in our expansion cannot endanger the stability of the 
strong-coupling fixed point. 

As can be deduced from Fig. 0, our approximations become less accurate for < Air and eventually our present 
approach breaks down for < 2tt (Sect. IV. C). In addition, the bound states present in the spectrum for < Att 
according to the exact solution are absent in our solution. These bound states will be generated by interactions in 
-ffios that are not included in our present approximation. Work on these issues for < Air is in progress. 



To summarize, we h ave obtained an explicit approximate relation between the strong-coupling problem (183) and 



its diagonalized form (184) without using the integrable structure of the model. The method presented here provides 
a theoretical tool that is capable of achieving this in a systematic expansion throughout the crossover region. We have 
been able to carry out this program completely for > An, and obtained first results for < An (where more work 
remains to be done e.g. regarding the bound states). Our approach is conceptually simple since a small parameter 
is identified and used as an expansion parameter. It is therefore possible to study nonintegrable perturbations 
with our approach, and the calculation of correlation functions also seems feasible. Finally, there are various other 
onedimensional strong-coupling problems, as for example the Kondo model, where the present approach should be 
useful ||. 

The author acknowledges many valuable discussions with W. Hofstetter and D. S. Fisher. This work was supported 
by the Deutsche Forschungsgemeinschaft (DFG), by the SFB 484 of the DFG and by the National Science Foundation 
(NSF) under grants DMR 9630064, DMR 9976621 and DMR 9981283. 

APPENDIX: PROPERTIES OF VERTEX OPERATORS 

1. Operator product expansion 

This Appendix compiles some important properties of vertex operators. For a review of these properties see also 
Ref. [|0). 

We first want to establish a relation between the normal-ordered and the non-normal-ordered vertex operators. 
By definition 

^\ -.-fbl-»PVi(p) I : 



Vi(a; X) — : exp a — e 

= exp U^^e-tlp|-»P-(ji(p) j cxp 

\ p>0 ^ / \ p<0 



p 



Mp) (A1) 



and next we can use the formula 

gAgB^gC/2gA+B (A2) 

with C = \A, B\ since C is a number and commutes with both A and B: 
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L p>o 

p>0 



CTl 



g<0 



-a|p| 



The sum over q yields 



and in the thermodynamic hmit L —f oo 



C = -a^ In 



27ra 

~17 



Therefore 



Vj- (a; a;j = : exp ±a 



\2TTa J 



exp 



p#0 



e 2 



An important property of the Fourier-transformed vertex operators ([78|) is their action on the vacuum 

Vi{a;-k)\n) ^ Vi{-a]k)\n) = V2{a]k)\VL) = V2{-a]-k)\VL) =0 Vfc > . 
This is shown easily, e.g. 



V^{-a-k)\n) = / dxe-'^"1/i(-a;x)|r!) 



1 

2^ 



dx e '''^ exp a ^ 

\ P>0 



One expands the second exponential and arrives at terms with the structure 

j dx e-"^''+P'+-+P'^^ = 

since . . . ,p„ > 0. Therefore 

Vi{-a;k)\n) =0 



for fc > and the analysis for the rest of (A7) proceeds likewise. 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 



(AlO) 



Next we want to want to evaluate expectation values for products of vertex operators. We look at the product 



Vi{a;x)Vi{-a;y) 



exp(-f IpI - ipx) 



exp(-f |p| - ipx) 



X exp j —a 



exp(-||g| - iqy) 



ai{q) j exp I 

/ \ 9<0 



exp(-f l^l - iqy) 



and commute the second and third exponentials using the formula 



A B c B A 
e e — e e e 



(All) 



(A12) 



with C = [A, B] since again C is a number and commutes with both A and B: 
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exp(-2 bl - tpx) ^ exp(-2 |g| - iqy) , , 
a > 7= fJib), ~a > j= CTi(g) 



P<0 V g>0 

2 exp(-a|g| + iq{x ~ y)) 



9>0 



In the thermodynamic hmit L oo this gives 



L/2tt 



i{y — x) + a 



Therefore 
Vi{a]x)Vi{-a;y) 



L/2tt 



i{y -x)+i 



2 / 



/ ^ cxp(-f|p| -ipx) \ / 
exp -a 2^ cri(p) exp 

V t^o v\p\ J \ 



(71 (p) 1 exp I 

/ \ g>o 

exp(-f |g| - iqy) 



cxp(-f Igl - iqy) 



'^1(9) 



(A13) 
(AM) 

(A15) 



p<0 

From this we obtain the expectation value 

ma;x)V,i-a;y)) 



L/2tt 



(A16) 



i{y — x) + a, 

A similar calculation can be done for V2{a;x) and the only difference is the exchange of x and y in the denominator 

L/2iT 



{V2{a-x)V2{-a-y)) 



i{x — y) + a 



(A17) 



The operator product expansion (OPE) for vertex operators can be deduced from ( A15 ) 

L/2tt 



Vi{a;x)Vi{-a;y) = 



i{y -x)+i 



-X expLE^^^^^(^-----)-iW 



exp j —a 

\ P<0 



exp(-fbl) 



(e- 



0-1 (P) 



L/27: 



i{y — x) + a 



exp -zaix -y)J2 e^^lPl-'^'/H ai(p) + 0((a; - y)^ 



p>0 



exp j —ia{x — y) ( 

V p<0 



e'^M-^P-^\ai{p)+0{{x-y)^) 



L/2tt 



i{y — x) + a 



1 + za(y - x) ^ e-tlPl-P-/^ai(p) + 0{{x - y)') 

p#0 



(A18) 



Higher order terms in the OPE can easily be deduced using the above scheme. These terms are less singular as a; ^ y, 
or, in the language of renormalization theory, they have a larger scaling dimension (are more irrelevant) and can be 
expressed as spatial derivatives of the bosonic field. 
A similar calculation for V2(a;x) gives 



V2{a;x)V2{-a;y) 



L/2tt 



i{x — y) + a 



1 + ^a{x -y)J2 e-^l^l-*^"/W ^2(p) + 0((x - y)' 

p#0 



(A19) 
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again the only difference is the exchange of x and y. 

It is well-known that for a = ±1 the vertex operators describe fermion creation and annihilation operators. This 
can be checked easily by using the OPE (A18) in the anticommutator for the special case a — 1 



{Vi{l;x),Vi{-l;y)} = 



L/2tt 



L/2tt 



i{y ~ x) + a i{x — y) + a 



1 + taiy -x)Y, e-tlPl-'^f^/H + 0((a; - yf) 

p#0 



L6{x-y) l + ia(y-:E)^e-tlp|-'P-/^ai(p) + 0((x-y)2) 
= L5{x~y) 

in the limit a ^ 0. All higher order terms in the OPE vanish in this limit. Likewise one finds 



{Vi{l;x),Vi{l;y)} = {Vi (-1; x), y)} 

and the same relations for V2(±l; x). 



(A20) 



(A21) 



2. Exchange relations 

Let us look at the commutation relation of vertex operators in momentum space. These are worked out in the following 
for general a. For simplicity we only consider Vi(±a;x), the calculation for V2(±a;a;) proceeds along similar lines. 
From the definition of the vertex operators one easily verifies the following relation 

V^{-a;x)V^ia;y) = t^^^£±^VAa;y)VA~a; x) . (A22) 

For large distances |a; — ?/| 3> a the coefficient can be well approximated by 

cos(7ra^) + isgn(x — y) sin(7ra^) , (A23) 

whereas for small distances it becomes equal to 1. Now for small distances the operator product expansion for the 
two vertex operators can be used and it is then possible to write generally for all x ^ y 

Vi{-a;x)Vi{a;y) - OPE(x ^ y) = (cos(7ra2) + isgn(j; - ?/) sin(7ra^)) {Vi{a;y)Vi{-a; x) - OPE(a; y)) (A24) 

For our purposes here it will be sufficient to look only at the leading c-number term in the OPE (higher orders can 
easily be taken into account if necessary). This is equivalent to subtracting the ground state expectation value on 
both sides (|50|). Notice that our relation "closes" 

* Vi (—a; x)Vi {a; y) * = (cos(7ra^) + zsgn(x — y) sin(7ra^)) * Vi {a; y)Vi (—a; x) * 
— (cos(7rQ;^) + isgn(a:: — y) sin(7rQ;^)) 

X (cos(7rQ!^) + isgn(y — x) sin(7ra^)) *Vi{—a;x)Vi{a\y)* 
= Wi{~a-x)Vi{a-y)* (A25) 

In terms of Fourier components (^8|) this reads {a > 0) 

*Vi{-a;ki)Vi{a;k2) * = -^'^*Vi{a; k2 + p)Vi{-a; ki + p) * f e*^^ (cos(7rQ:^) + z sgn(x) sin(7ra^)) 

= cos(7ra^) * Vi{a\ k2)Vi{—a; ki) * 

2 °° 1 r / 27r \ / 27r \ 

+ sin(7ra2) -V *Vi [ a;k2 - —{2n + 1)] Vi -a; fci - — (2n + 1) * 



Vi I a;A:2 + ^(2n + l) ) Vi ( -a; fci + ^(2n + 1) ) * 
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cos(7ra^) * Vi{a; k2)Vi{—a; fci) * 

1 f 1 

sin(7rQ!^) Re / dq * Vi{a; k2 + q)Vi(—a; ki + q) * 

TT J q + le 



TT 



dq 



1 



q + ie 

in the thermodynamic hmit with hme^o being understood. Likewise 



* Vi{a; k2 + q)Vi{-a; fci + q)* 



* Vi{a; ki)Vi{~a; k2)* 



1 



Im 



dq ■ 



1 



q + ie 



Vi{-a; k2 + q)Vi{a\ fci + q)* 



(A26) 



(A27) 



Notice that these relations close in fc-space in the same sense as ( A25 ). They are only "simple" for integer when 
the vertex operators behave as bosons or fermions: Then the term proportional to sin(7rQ;^) with its summation over 
wavevectors vanishes. 
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